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Abstract. - We study the effects of a high-frequency (HP) signal on the response of a noisy 
bistable system to a low- frequency subthreshold sinusoidal signal. We show that, by conveniently 
choosing the ratio of the amplitude of the HF signal to its frequency, stochastic resonance gains 
greater than unity can be measured at the low-frequency value. Thus, the addition of the HF 
signal can entail an improvement in the detection of weak monochromatic signals. The results are 
explained in terms of an effective model and illustrated by means of numerical simulations. 



The phenomenon of stochastic resonance (SR) has 
been studied with growing interest during the last three 
decades, being found to be of relevance in a great variety of 
phenomena in physics, chemistry, and the life sciences [1]. 
Roughly speaking, SR consists in the amplification of a 
weak, time-dependent signal of interest by the concerted 
actions of noise and the nonlinearity of the system. Sev- 
eral quantifiers have been used to characterise the SR re- 
sponse of noisy systems in the presence of periodic signals. 
In particular, the nonmonotonic behaviour of the output 
signal-to-noise ratio (SNR) with the strength of the noise 
is a widely accepted signature of SR. In addition, a dimen- 
sionless quantity known as the SR gain is usually defined 
as the ratio of the output SNR over the input SNR. The 
SNR measures the "quality" of the signal, in terms of the 
ratio of its "coherent" (periodic) component over its "inco- 
herent" (noisy) component. In turn, the SR gain compares 
the "qualities" of the output and the input signals. In gen- 
eral, obtaining high output SNRs and SR gains greater 
than unity would be desirable when using SR as an ampli- 
fication mechanism. Analog [2,3] and numerical [4,5] sim- 
ulations have shown that noisy bistable systems driven by 
subthreshold multifrcqucncy forces can display SR gains 
greater than unity when the parameters of the problem 
are properly chosen. Moreover, in [6], a two-state model 
of SR has been used to explain these results analytically. 
By contrast, to the best of our knowledge, there is no evi- 
dence of SR gains greater than unity when a subthreshold 
monochromatic signal drives an isolated bistable system. 



However, these unusual large gains have been reported in 
the case of a suprathreshold sinusoidal signal for an iso- 
lated bistable system [7] and in the case of a subthreshold 
sinusoidal signal for coupled bistable systems [8] . 

As it has been already pointed out in the literature [7] , 
an improvement of the response of a nonlinear system to 
a subthreshold signal of interest embedded in noise can 
be achieved better by lowering the threshold value rather 
than by increasing the noise strength. However, in most 
cases of practical interest, threshold lowering is usually a 
difficult task. Recently, it has been shown that the ef- 
fect of a strong high-frequency (HF) monochromatic force 
on the overdamped dynamics of a Brownian particle in a 
bistable potential can be described in terms of an effective 
potential whose characteristics depend on the parameters 
of the HF field [9]. In particular, the barrier height of this 
effective potential is smaller than the original one. Thus, 
an HF force provides a mechanism to achieve an effective 
threshold lowering. As shown in [10] for a square-wave 
signal of interest, it is possible to take advantage of this 
effective threshold lowering to improve the SR gain values 
evaluated at the fundamental frequency of the square wave 
signal. More precisely, in [10], SR gains greater than unity 
are found for a wide range of values of the noise strength. 
However, as mentioned in the same reference, this positive 
effect of the HF force seems to be absent when the signal 
of interest is sinusoidal. It is pertinent to point out that 
the effect of HF fields on nonlinear stochastic systems is 
not necessarily positive (see, for instance, [11] for a ratchet 
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model and [12] for an excitable model). Nevertheless, in 
ratchet models positive effects have been also found [13]. 

In this work we show that, by using a strong HF sig- 
nal, SR gains greater than unity can be also observed in 
a bistable noisy system when the signal of interest is sub- 
threshold and sinusoidal. In order to do that, let us con- 
sider a stochastic system characterised by a single degree 
of freedom x{t), whose dynamics (in dimensionless units) 
is described by the stochastic differential equation 



x{t) = -U'[x{t)]+F{t)+^{t). 



(1) 



Herein, ^{t) is a Gaussian white noise of zero mean and 
autocorrelation function {^{t)^{s)) = 2DS{t— s), D being 
the noise strength, F{t) is the sinusoidal driving F{t) = 
A cos(r2t), and U'{x) denotes the derivative with respect to 
X of the quartic potential U{x) = a;^/4— x^/2. Henceforth, 
we will restrict ourselves to subthreshold driving signals. 
More precisely, we will assume that A < Ath = 2/-\/27, 
where ^th is the static threshold value. In this case, the 
time dependent potential U{x) — Ax cos{i}t) possesses two 
local minima separated by a local maximum for any in- 
stant of time t. Notice that the dynamical threshold value 
(defined as the max;imum value of A such that no inter- 
well transitions are possible in the absence of noise) always 
exceeds this adiabatic threshold ^th- 

Equation (1) describes the overdamped motion of a par- 
ticle in a symmetric double well potential driven by a sinu- 
soidal force and noise. Alternatively, x('")(t) = F{t)+^{t) 
may be interpreted as an input signal consisting of a noisy 
term ^{t) and a periodic term F{t) (signal of interest). 
This input signal is transformed into an output signal 
a;(°"t)(t) = x{t) after being processed by the nonlinear 
device characterised by U{x). In this context, SR can be 
understood as a mechanism in which noise plays a positive 
role in the optimal detection of some features of the signal 
of interest F{t) in the output signal x{t). 

Several quantifiers have been used to characterise the 
SR phenomenon. Here, we will only consider those quan- 
tifiers involving the power spectral density (PSD) of the 
input and output signals, namely, the SNR of the output 
signal and the SR gain. According to the generalisation 
of the Wiener-Khinchine theorem to a periodically driven 
stochastic process, the PSD of the output signal x{t) is 
given by [14] 



2 r°° 

5(°"*)(a;) = - / drC( 



r) cos(a'T), 



(2) 



with Lo > 0, where C(r) = dt{x{t + T)x{t))oo, 

{' ' ■)oo representing the average over the realizations of the 
noise after a relaxation transient stage. It can be proved 
that S'^™*) (u>) consists of a scries of deltalike spikes located 
at the odd harmonics of the fundamental frequency fl, 
superimposed on a background PSD [1]. This background 
PSD (also called incoherent part of the PSD) is given by 
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where Cincoh(T) = C(t) - Ccoh(T), with Ccohir) = 
Jq dt {x{t + T))oo{x{t))oD (see, for instance, [15]). 
Then, the output SNR is defined as [1,15] 
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The input PSD 5(^")(w), as well as slZl^{u;) and 
can be obtained from the above expressions by replacing 
x{t) by a;(™^(t). In particular, it is easy to prove that the 
input SNR reads 



AD 

Finally, the SR gain is defined as 



G 



^(out) 



(5) 



(6) 



As mentioned before, obtaining SR gains greater than 
unity would be desirable. For this purpose, let us intro- 
duce a strong HF signal of the form: 



Y{t) = NQr cos{Nnt + if) , 



(7) 



where the parameter N is a. positive integer, r is the ratio 
of the amplitude of Y{t) to its frequency, and is an 
arbitrary initial phase. We are interested in situations 
in which the parameters Nflr and Nfl appearing in this 
monochromatic force are much larger than the rest of the 
parameters in the problem. That is why Y(t) is called a 
strong HF signal. This situation can be formally achieved 
by taking the limit N oo, with the ratio r kept fixed. 
The HF signal Y{t) must be summed to the right hand side 
of eq. (1), so that now x{t) fuUfils the Langevin equation 



x{t) = -U'[x{t)] + F{t) + Y{t) + e(t) . 



(8) 



As a consequence, the input signal in our previous discus- 
sion must be replaced by x^'^'Xt) = F{t) + Y{t) + ^{t). 
It is easy to prove that the input PSD in the presence of 
the HF signal [5'('")(w;r) with r ^ 0] can be expressed in 
terms of that obtained in its absence [S'('")(a;; 0)] as 

(w; r) = (w; 0) + d{LJ - Nfl) . (9) 

It should be emphasised that we are interested in the eval- 
uation of the input and output SNRs at the frequency fl 
[frequency of the signal of interest F{t)], Y{t) being just 
a tool introduced to improve the response features of the 
system at that frequency. In addition, since in the defini- 
tion of only the PSD around the frequency fl plays 
a role, it is clear from eq. (9) that the input SNR is given 
by eq. (5) even in the presence of the HF signal Y{t), as 
in this case iV 1. 

For reasons that will be clarified later, it is convenient to 
define the new stochastic process z{t) = x{t) — r sm{Nrtt + 
ip). After replacing this definition in eq. (2), it can be 
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proved that the output PSD, and the PSD of 

the process z{t), Sz{co), are related according to 



100 



5(°"*)(w) = 5^(w)+2r 



-lm(e-'^Zjv) S{(j-Nn), 
(10) 

where Zn is the Fourier component of the process z{t) at 
the frequency Nil, i.e., Zn = T'^ dt{z(t))ooe-^^^K 
From cq. (10) it is clear that, in order to evaluate quan- 
tities which only depend on the behaviour of the output 
PSD around the frequency [as for instance both 
processes, x{t) and z(t), are equivalent as long as N ^ 1. 

The dynamics of the stochastic process z{t) is described 
by the equation 

z{t) = -U'[z{t) +rsm{Nnt + if)] + F{t) + ^{t) . (11) 

As we have mentioned before, we are interested in the 

limit N oo, with r kept fixed. From eq. (11), it follows 
that the time derivative of the process z{t) is at most of 
order 1 as A'' ^ oo. In this sense, we will say that this 
process is slow. In contrast, according to the definition 
of z{t), it is clear that, in this limit, the process x{t) is 
fast, in the sense that it is a highly oscillating process 
around z{t). As the process z{t) is slow, a large number 
of oscillations of the function r sin( Afii + (p) appearing in 
U'[z{t) + rsm{NO,t + (p)] takes place before a significant 
change in z(t) occurs. Therefore, for A ^ 1, z{t) must 
be almost independent of the phase (p and, consequently, 
z{t,(p) ~ z{t) := {27r)~^ f^^ dip' z{t, if'), where we have 
written explicitly the dependence of z{t) with respect to 
(p. Taking into account this property and carrying out the 
phase average in eq. (11), one obtains 



where we have introduced the effective potential 
1 /-^^ 

Ues[z] ■■= — / dpU[z + rsm{Nnt + p)] 

271" Jo 



(12) 



—4 -2 2 



3r2 



- 1 



(13) 



with c(r) = 1 — 3r^/2. From this result, it is clear that 
the stability of this effective potential depends on the ratio 
r. More precisely, if r < s/^J^ the potential is bistable, 
whereas if r > \f^j2> is monostable. Thus, an increase in 
r leads to a decrease in the effective barrier height and, 
eventually, to its disappearance. 

Henceforth, we will restrict our study to the case r < 
y^2/3, which corresponds to < c(r) < 1 (bistable poten- 
tial). In order to reduce the effective Langevin equation 
[eq. (12)] to the standard form [eq. (1)], it is convenient to 
introduce the following time and space scale changes: 



i = tc{r), 
z{i) = z[i/c{r)]/^/^ 



(14) 
(15) 



10 





1 




^2 




fJR^,£l)=[R^,C}/c(r)] 








^.^y-^;\R^,ac(r)]=(R^,a) ; 




R 

1 

. ...1 


1 1 . . 





101 



0.01 



0.1 
D 



10 



Fig. 1: Two regions, Ri and R2, of the parameter plane 
D — A are transformed into each other by the mappings fr and 
f~^. Ri is the subthreshold region delimited by the four ver- 
tices (75, A/Ah) = (0.005,0.1), (0.015,0.1), (0.015,0.3), and 
(0.005,0.3), and R2 is obtained from Ri by applying fr with 
r = 0.800529. This last parameter value has been chosen so 
that R2 lies within a region wherein SR gains greater than 
unity have been observed in the case of a bistable system in the 
presence of a suprathreshold sinusoidal signal with frequency 
n = 0.1 (see fig. 2 in [7]). 



Then, it is easy to prove that the stochastic process z{i) 
fuUfils the Langevin equation (1) but for the rescaled val- 
ues of the amplitude, noise strength, and frequency given 
by ^/c(r)'^/^, _D/c(r)^, and 0/c(r), respectively. We will 
denote by the mapping which assigns to each set of pa- 
rameter values {A, D, f2) the corresponding rescaled val- 
ues, and by f~^ its inverse mapping, i.e., 



fr{A,D,n) 

f-\A,DM) 



A/c{rf'\Dlc{rf,n/c{r)\,{m) 
Ac(rfl^,Dc{rf,nc{r)\. (17) 



In fig. 1 we depict how two regions of the parameter plane 
D — A, Ri and R2, are transformed into each other by the 
mappings fr and f~^. 

As mentioned before, for N ^ 1, the output SNR of 
the processes x{t) and z{t) coincide exactly. Furthermore, 
for A » 1, the processes z{t) and z{t) give rise to the 
same value of this quantifier with a great deal of approxi- 
mation. Finally, the process z{t) is related to the process 
z{i) [which is the process x{t) but in the absence of the 
strong HF signal and for the above mentioned rescaled pa- 
rameter values] by means of the transformations (14) and 
(15). Then, by using the definition of the output SNR, 
one obtains that 

[{A, D, O); r] ~ c(r)i?(°"*) [/.(A D, fi); 0] (18) 

for A » 1, as expected from the change of scale in eq. (14) 

and the fact that the SNR has dimensions of timc"-'^. The 
above expression relates the output SNR in the presence of 
the HF signal, i?(°"*)[(^,£),n);r], with that correspond- 
ing to a new set of parameter values but in the absence of 
the HF signal, [fr{A, D, Vl); 0]. 
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Fig. 2: Dependence on the noise strength D of the SR gain G 
in a bistable noisy system driven by a sinusoidal subthreshold 
signal, both in the absence (empty diamonds) and presence 
(full diamonds) of an HF signal. The parameter values of the 
sinusoidal subthreshold signal F{t) and the HF signal Y{t) are: 
A = 0.1, n = 0.005, r = 0.797, N = 4000, and (p = 0. The 
SR gain obtained from the numerical solution of the effective 
dynamics in eq. (12) is plotted with solid line. The horizontal 
dashed line marks the unity for the SR gain. 



Analogously, taking into account the mapping in 
eq. (16), as well as the fact that the input SNR at the 
frequency of interest Q is given by eq. (5) even in the pres- 
ence of the HF signal [see the paragraph below eq. (9)], it 
is also clear that 

i?(-)[(AA^);r] = R^"'^[{A,D,n);0] 

= c(r)i?('") [/,. (A, i:>,f7);0] .(19) 

Finally, from eqs. (18), (19), and the definition (6), one 
obtains that 



G[iA,D,n):r]c^G[fr{A,D,n);0] 



(20) 



for iV 1, this last equation being just a consequence of 
the dimensionless character of the SR gain. Equation (20) 
is one of the main results of this paper. It shows that the 
value of the SR gain in the presence of an HF signal is 
approximately equal to that corresponding to a system in 
which that HF signal is absent, but for a different set of 
parameter values. Obviously, the larger the value of N, 
the better the approximation. 

From eq. (20) it is easy to prove that, by applying 
an appropriate strong HF signal to a bistable noisy sys- 
tem driven by a subthreshold sinusoidal signal of inter- 
est, one can obtain SR gains greater than unity. Indeed, 
as shown in [7], SR gains greater than unity can be ob- 
served in a noisy bistable system driven by a suprathresh- 
old sinusoidal signal. For instance, in fig. 2 of [7] it is 
shown that, for a fixed value of the frequency of this sig- 
nal O = 0.1, there exists a region of the parameter plane 
D — A wherein G > 1. To fix ideas, let us assume that 
G[{A2, D2, O2); 0] > 1, with A2 > Ath, and let us choose a 



value ro such that J2 [l - {Ath/A2)^/^] /3 < ro < ^/2/3. 



If we define {Ai, Di,Cli) = f^^^{A2, 02,^^2), it is easy to 
check that Ai < Ath- Furthemorc, according to eq. (20), 
G[(Ai,£)i,l^i);ro] ^ G[(A2, i'a, 1^2); 0] > 1, for iV » 1. 
Thus, by applying an HF signal of the type given in eq. (7) 
(with f2 = lli, r = ro, iV^ 1, and an arbitrary value of 
(f) to a noisy bistable system driven by a subthreshold 
sinusoidal signal of frequency fli and amplitude Ai, one 
gets an SR gain greater than unity at a noise strength Di. 
This is exactly what we wanted to prove. 

As an example of the above reasoning, let us consider 
fig. 1. There, R2 lies within a region of the parameter 
plane D—A wherein SR gains greater than unity have been 
observed in the case of a bistable system in the presence of 
a suprathreshold sinusoidal signal with frequency Q = 0.1 
(see fig. 2 in [7]). Therefore, according to our reasoning, it 
is expected to observe SR gains greater than unity within 
the subthreshold region i?i and for the frequency value 
Q = 0.003873, whenever one applies an HF signal of the 
type given in eq. (7), with r = 0.800529, A/' > 1, and an 
arbitrary value of (p. 

In order to illustrate the feasibility of the above argu- 
ment, in fig. 2 we depict the dependence on the noise 
strength D of the SR gain in a bistable noisy system driven 
by a sinusoidal subthreshold signal, both in the absence 
(empty diamonds) and presence (full diamonds) of an HF 
signal. The SR gain has been evaluated numerically fol- 
lowing the method described in [4]. The parameter val- 
ues of the sinusoidal subthreshold signal are A = OA and 
ri = 0.005. As expected, in the absence of the HF signal, 
the SR gain is lower than unity for all the noise strength 
values. Nevertheless, when an HF signal of the type given 
in eq. (7) is applied, with n = 0.005, r = 0.797, N = 4000, 
and (fi = 0, there is a range of noise strength values for 
which SR gains greater than unity are observed, the max- 
imum value of G being located aroimd -Dmax = 0.01. Since 
our argument is premised on the equivalence of the exact 
and the effective dynamics [eqs. (8) and (12), respectively] 
for evaluating the SR quantifiers in the asymptotic limit 
A'' ^ 00, it is important to check the validity of this ap- 
proximation for N = 4000. In fig. 2, we have also plotted 
with solid line the SR gain obtained from the numerical 
solution of the effective dynamics in eq. (12). A glance 
at that figure reveals that the agreement between both 
dynamics is excellent for A^ = 4000. 

It is worth noting that the rescaled parameter values 
of the monochromatic signal in fig. 2 [i.e., A/c{r)^^'^ « 
25.347^th and i}/c{r) ss 0.106] are almost equal to the 
amplitude and frequency of the suprathreshold sinusoidal 
signal of the dot-dashed line in fig. 5 of [7]. Thus, accord- 
ing to eq. (20) ; the above mentioned solid and dot-dashed 
lines must be essentially the same, except for a rescaling of 
the D-axis by a factor l/c(r)^ « 449.122. So, for instance, 
the location of the maximum of the SR gain in the dot- 
dashed line is approximately equal to r>max/c(r)-^ ~ 4.491, 
and the maximum value of the SR gain is approximately 
1.2 in both curves. 

In conclusion, in this work we have shown that when 
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a noisy input containing a subthreshold monochromatic 
signal of interest is processed by a bistable system in the 
presence of an HF signal, SR gains greater than unity 
can be observed. In order to explain this result, we have 
proved that the original bistable system in the presence 
of the HF signal can be approximated by a new bistable 
system obeying the same stochastic dynamics as the orig- 
inal one but in the absence of the HF signal. In this new 
system, the amplitude and frequency of the signal of in- 
terest, as well as the noise strength, appear rescaled by 
coefficients depending on the ratio r of the amplitude of 
the HF signal to its frequency. The main result of the 
paper establishes that the SR gains of both systems are 
approximately equal. The importance of this result relies 
on the fact that, although the original parameter values 
lie within a subthreshold region, the rescaled parameter 
values obtained from a convenient choice of r can belong 
to a suprathreshold region wherein SR gains greater than 
unity have been observed. As the SR gains in both sys- 
tems approximately coincide, it is apparent that the ad- 
dition of the HF signal can give rise to SR gain values 
greater than unity. Finally, the theoretical results have 
been illustrated in a particular case by means of numeri- 
cal simulations. The authors are confident that the results 
reported in this paper can be useful for the optimal detec- 
tion of weak monochromatic signals in the great variety of 
systems modeled by noisy bistable potentials. 
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